Is the particle current a relevant feature in driven lattice gases? 
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By performing extensive MonteCarlo simulations we show that the infinitely fast driven lattice gas 
(IDLG) shares its critical properties with the randomly driven lattice gas (RDLG) . All the measured 
exponents, scaling functions and amplitudes are the same in both cases. This strongly supports the 
idea that the main relevant non-equilibrium effect in driven lattice gases is the anisotropy (present 
in both IDLG and RDLG) and not the particle current (present only in the IDLG). This result, 
at odds with the predictions from the standard theory for the IDLG, supports a recently proposed 
alternative theory. The case of finite driving fields is also briefiy discussed. 

PACS numbers: 64.60.-i, 05.70.Fh 



The Ising model exhibits a prototypical equilibrium 
phase transition, and the associated cj)'^ Ginzburg-Landau 
theory is a paradigm of continuous theory for equilib- 
rium critical phenomena However, thermodynamic 
equilibrium is exceptional in nature where stationary 
states are typically away from equilibrium With 
the purpose of defining simple lattice models describing 
generic non- equilibrium phase transitions, different ten- 
tatives have been made in the last two decades. Among 
them, perhaps the most intriguing example is the driven 
lattice gas (DLG) |^-^. (Other interesting examples 
are the directed percolation model j|] and the Kardar- 
Parisi-Zhang equation). The DLG, being a straightfor- 
ward extension of the Ising model, has in fact become 
a workbench for emergent non-equilibrium theories and 
field theoretical approaches. 

The DLG is a d-dimensional kinetic Ising model with 
conserved dynamics, in which transitions in the direction 
(against the direction) of an externally applied field, E, 
are favored (unfavored) while transitions perpen- 

dicular to the field are unaffected by it. The field in- 
duces two main non-equilibrium effects: (i) the presence 
of a net current of particles along its direction, and (ii) 
strongly anisotropic configurations . At high tempera- 
tures, the system is in a disordered phase while, for half- 
filled lattices (the only case we refer to in what follows) 
there is a second-order critical point, below which the 
DLG segregates into (two) high and low density aligned- 
with-the-field stripes. Establishing unambiguously the 
DLG universality class is an important issue in the way 
to rationalize the behavior of non-equilibrium systems. 

Continuous approaches such as Langevin and associ- 
ated field theories [|| have been most useful in study- 
ing universality issues in equilibrium critical phenomena. 
In particular, coarse-grained approaches combined with 
renormalization group (RG) techniques provide a method 
for the classification of the different possible terms (op- 
erators) as relevant, irrelevant or marginal. In fact, 
Langevin equations are more illuminating than other 
(even more rigorous) approaches, as they permit to un- 



derstand systematically how possible perturbations or 
model variations would affect critical properties. Con- 
sequently, many studies have focused on the DLG and 
its universality by using both non-equilibrium continu- 
ous approaches and computer simulations (unfortunately, 
general exact solutions are not available). Within this 
perspective, it is somewhat deceptive that after many 
computer and analytical studies, the universality class of 
the DLG remains a debated issue ■ 

A phenomenological Langevin equation intended to 
capture the relevant physics of the DLG at criticality 
was proposed and renormalized more than a decade ago 
Q. This equation, referred to as driven diffusive sys- 
tem (DDS), is a natural extension of the conserved cj)'^ 
theory for the Ising equilibrium transition (model B Q|) 
and seems to capture the main symmetries and conserva- 
tion laws of the discrete DLG. It includes a particle cur- 
rent term (which from naive power counting turns out to 
be the most relevant nonlinearity) as well as anisotropic 
coefficients. It certainly is a suitable and very reason- 
able candidate to be the canonical coarse-grained model 
representative of the DLG universality class. The DDS 
Langevin equation reads: 
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where (j) is the coarse grained field, C is a conserved Gaus- 
sian noise and the cubic term, being a dangerously irrel- 
evant variable is kept in order to ensure stability 
T||, T±, A and a are model parameters. The most em- 
blematic (exact) prediction derived from the DDS RG- 
analysis, namely, the mean field behavior of the order pa- 
rameter critical exponent, /3 = 1/2, ||] has eluded a large 
number of Monte Carlo (MC) analysis aimed at prob- 
ing it jl2|, however. In particular, systematic deviations 
from the predicted scaling are observed both in d = 2 
|,ra and in d = 3 Indeed, different MC analy- 

sis (performed using a variety of aspect-ratios and order 
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parameters) lead systematically to a value of /? close to 
~ 0.3 (in d = 2), with error bars excluding apparently 
the value /? = 1/2 (see ||] for a critical review of sim- 
ulation analysis). This is a main indication that, strik- 
ingly enough the DDS equation does not describe properly 
the infinitely fast driven DLG (IDLG) critical properties. 
Moreover, there are some other hints suggesting strongly 
that the differences between the predictions of the stan- 
dard Langevin approach and MC results are more fun- 
damental than a simple discrepancy in the value of (3. In 
particular, the intuition developed from MC simulations 
of the DLG and variants of it ||^ suggests that, contrarily 
to what the DDS equation establishes, it is the anisotropy 
and not the presence of a current the basic ingredient con- 
trolling the critical behavior For instance, in a modi- 
fied DLG in which anisotropy is included by means others 
than a current , the scaling behavior at criticality re- 
mains unaltered upon switching on an (infinite) driving 
(see |l3|,^). Other compelling evidences supporting this 
hypothesis can be found in P,p^. 

In an attempt to clarify this puzzling situation, and 
reconcile continuous approaches with numerics, different 
scenarios have been explored. In particular, an alterna- 
tive route to build up Langevin equations starting from 
generic microscopic master equations was recently pro- 
posed . By applying this approach to the DLG, one ob- 
serves that, owing to a transition-rates saturation effect, 
the coefficient a of the non-linear current term, V||(/)^, 
vanishes in the limit of infinite driving fields and, there- 
fore, it does not appear in the final Langevin equation nor 
it is generated perturbatively [Q. The resulting theory 
(alternative to Eq.(P) is: 

dtCj^{v, t) = T^^lcj, - S/^cj, + ^Vi03 + T|| Vjj0 + C (2) 

plus higher order irrelevant contributions (note that a 
linear current term has been eliminated by employing a 
Galilean transformation [^,0). This equation, named be- 
low anisotropic diffusive system (ADS), is a well known 
one: it coincides with the Langevin equation representing 
the random DLG (RDLG) ||l5| , [l6|] (for which the driving 
field takes values oo and — oo in a random unbiased fash- 
ion, generating anisotropy but not an overall current). 
This theory has been extensively studied in [ p^p6| ; its 
critical dimension is = 3 (instead d^ = 5 for the 
DDS) and the critical exponents and finite size scaling 
(FSS) properties are now well known. Other systems 
in this universality class are the two-temperature model 
and the ALGA model This theory for the IDLG 
includes anisotropy as its basic non-equilibrium ingredi- 
ent. Instead - for non-saturating, finite, driving fields - 
the cancellation of the nonlinear current term does not 
occur, and our method recovers the standard DDS equa- 
tion. 

Aiming at further clarifying these issues, we report 
here on extensive MC simulations of the IDLG and the 
RDLG in d = 2. The main objectives are: (i) trying to 
conclude whether the IDLG and the RDLG share the 



same critical behavior or not; and (ii) measuring the 
critical exponents by performing systematic anisotropic 
finite-size scaling (FSS). In fact, we perform FSS analy- 
sis for both the IDLG and the RDLG by following the 
anisotropic FSS scheme proposed in | [T5[ | consistent with 
the ADS theory; this allows us to analyze systematically 
possible scaling differences between both models. We also 
report on the case of finite-driving DLG. 

We consider rectangular lattices of size L|| x L± with 
periodic boundary conditions and random sequential up- 
dating ^,||; the external field E acts in the x (parallel) 
direction. Particles jump to a randomly chosen nearest 
neighbor site (provided that it is empty) with probabil- 
ity: min(l, exp[— /3(Ai7 -I- EAj)]), where AH is the en- 
ergy (Ising Hamiltonian) variation, and Aj = (—1,0,1) 
for jumps along, against, and orthogonal to the direction 
of the field, respectively. Following ||l^,0 the order pa- 
rameter is chosen as the structure factor S'(0, 27t/L±). In 
order to perform a systematic anisotropic FSS we consid- 
ered system sizes 20 x 20, 45 x 30, 80 x 40, and 125 x 50. 

These aspect ratios satisfy Ljj^^'^" = 0.2236 x L±, where 
iy^/zy|l « 1/2 consistent with ADS anisotropic spatial 
scaling ||l^,|l^ . The number of MC steps considered var- 
ied between 1.8 x 10* and 2.4 x 10*, much larger than 
in any previously reported MC simulations. The total 
CPU time employed is about eight months in a pentiu- 
mlll 400Mhz machine. The critical temperature is de- 
termined by using the fourth (Binder) cumulant method 
JlSf . For the IDLG, the critical temperature is found 
to be Tf = 1.396(4)To (To is the Onsager temperature), 
slightly below previously reported values |^,^, while we 
find = 1.390(4)ro for the RDLG (see insets of Fig.3). 
These critical values were employed for the FSS analy- 
sis. In Fig.l we plot the order parameter, rescaled by 

a factor -^j^^''" , versus , where e is the distance to 

the critical point, for different system sizes A nearly 
perfect data-collapse is obtained by fixing = 1.25 and 
(3 — 0.33. The collapse gets worse upon slightly chang- 
ing these values; more precise estimates of the associated 
error-bars is a difficult and not essential issue in this con- 
text. Notice that we are plotting in the same graph data 
for the IDLG and for the RDLG, implying that the FSS 
scaling function is precisely the same for both models. 
Furthermore, the slopes of the asymptotic branches are 
approximately 1/3 and —0.61, consistent both with the 
order parameter exponent being f3 « 0.33, and 7 w 1.22 
(see below) [|l9|. In general, even when the scaling func- 
tions are universal their corresponding amplitudes are 
not expected to be so. For this reason, usually one has 
to introduce the, so called, metric factors (varying am- 
plitudes) [po[ in order to obtain superposition of scaling 
functions within the same universality class. Contrar- 
ily to this expectation, the magnetization scaling func- 
tions for IDLG and RDLG overlap perfectly. Therefore, 
it comes as a surprise that not only the scaling functions 
and the (3 exponent are universal in both models, but 
even the amplitudes coincide. A similar situation has 
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been recently reported for a different type of anisotropic 
FSS 111. 
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FIG. 1. Log- log plot of the order parameter rescaled by 
ij^^ " vs. eh^l^ " , for different system sizes: 045 x 30, 
A80 X 40, 0125 X 50. Filled (empty) symbols stand for the 
RDLG (IDLG); errorbars are smaller than the symbols. 
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FIG. 2. Log-log plot of the susceptibility rescaled by 
L|l vs. eL|| . Symbols are as in Fig.l (larger than er- 
rorbars) . 
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We have also computed the system susceptibilities, de- 
fined as the relative fluctuations of the order parameter: 
' < IT? > ~ < ra >^]. In Fig. 2, we plot the 

as a function of the rescaled 

II 



susceptibility times L|| 

distance to the critical point, ■ The best data col- 

lapse is obtained by employing the values 7 = 1.22 and 
i^ll — 1.25 for both models with, again, coinciding ampli- 
tudes. It should be stressed that this is the first time a 
really good collapse is observed below the critical point 
for anisotropic scaling of the IDLG. Plotting the dimen- 
sionless Binder cumulant as a function of the rescaled dis- 
tance to the critical point with = 1.25, again, nearly 
perfect data collapse is obtained for both models and all 
system sizes (Fig. 3). We also performed simulations in 
square lattices (128 x 128) as in some previous studies 
[^,|^. Monitoring m^^^ as a function of T/T^ we see 
no appreciable systematic difference between the curves 
for IDLG and RDLG, that have the same slope within 



numerical accuracy. The best linear fit correlation is ob- 
tained for P « 0.33 in both cases, providing an extra con- 
sistency check for our results. Moreover, eye inspection of 
IDLG and RDLG configurations, for any geometry, at a 
fixed relative temperature, does not permit to distinguish 
one from the other. In particular, the interfacial proper- 
ties look identical. Let us also stress that all the obtained 
exponent values are compatible with previous measures 
for the RDLG, as well as with the exponents obtained 
within an e— expansion of the ADS theory |l|jl|]. 




FIG. 3. Scaling plot of the fourth cumulant vs. eL^^^ " . 
Upper (lower) inset: fourth cumulant for the IDLG (RDLG) 
vs. T/To- Symbols are as in Fig.l. 

In conclusion, MC results support strongly that both 
the IDLG and the RDLG belong in the same universality 
class, and share not only critical exponents and scaling 
functions, but also the scaling amplitudes. This universal- 
ity class is described by the ADS equation, Eq.(||). There 
is absolutely no hint of any difference in the asymptotic 
behavior between the model with a current (IDLG) with 
respect to the current-less one (RDLG). All the numeri- 
cal evidence confirms that it is the anisotropy and not the 
net current (if any) the most relevant non- equilibrium in- 
gredient of driven systems. As discussed in the introduc- 
tion, this is striking from a field theoretical perspective 
given that the nonlinear current term, Vy^^, is naively a 
relevant perturbation at the ADS fixed point. In an alter- 
native approach, the coefhcient of such a term vanishes. 
In this picture, the fast drive limit corresponds to a sort 
of multicritical point in which an a priori relevant oper- 
ator is absent due to a cancellation of its coefficient and, 
consequently the usual "up-down" Ising symmetry (i.e. 
the three-point correlation functions vanish) is restored 
at criticality. In any case, it should be stressed that, from 
a more general perspective, field theoretical descriptions 
of non-equilibrium systems are much more delicate and 
subtle than their equilibrium counterparts, and an ex- 
tremely careful inspection of the system symmetries, con- 
servation laws, and dynamical features is required before 
venturing to make predictions on noncquilibrium univer- 
sality issues. For example, fermionic and bosonic non- 
equilibrium systems with the same dynamics, symmetries 
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and conservation laws have recently reported to belong 
to different universality classes [^. One could wonder 
whether the DLG hard-core interaction should be taken 
into account into its Langevin representation. 

Elucidating the critical behavior for finite E remains 
a challenging and interesting objective. Both, our al- 
ternative Langevin-building approach and the standard 
one, Eq.(|I|), include a relevant current term in this case 
and, consequently, predict [3 = f/2. Obtaining clear-cut 
results in this case is a computationally expensive task 
since: (i) As the external field appears in the argument 
of an exponential, even relatively small values of E gen- 
erate situations close to saturation, and strong crossover 
effects could hide the true asymptotic regime, (ii) If the 
field value is taken too small, crossovers from the equi- 
librium regime may also burden observations. A possible 
scenario that could follow from MC analysis is that finite 
fields show mean field behavior; that would be a strong 
backing for our theory Q that predicts the finite and 
the infinite driving cases to be qualitatively different. If, 
instead, scaling happens to be that of the ADS (as our 
preliminary MC results for i? = 3 and E — 0.5 seem to 
indicate; for E — 0.25 results do not quite fit this indi- 
cation) it would prove that it is for any arbitrary value 
of the driving field that anisotropy is the most relevant 
ingredient of driven systems. This scenario would un- 
cover a new puzzling situation and would certainly call 
for deeper theoretical understanding. Huge and careful 
simulations would be required to extract neat conclusions 
overcoming difficulties (i) and (ii) above. 

Summing up, we have performed extensive MC simula- 
tions of the IDLG and the RDLG. By using anisotropic 
finite size scaling techniques we have shown that both 
models belong to the same universality class: their crit- 
ical exponents, scaling functions and amplitudes are 
undistinguishable and coincide with those of the ADS 
equation. This result supports the conclusion that it is 
the presence of anisotropic coefficients, and not the par- 
ticle current the most relevant ingredient in these non- 
equilibrium driven problems (at least in the fast drive 
limit). Further theoretical efforts are certainly required 
in order to (i) to sort out if our alternative Langevin ap- 
proach is correct and what are its possible limitations, 
and (ii) to further clarify the universality issues of this 
quintessential non-equilibrium problem. Finally, it would 
also be very interesting to combine the powerful finite size 
methods recently introduced in this context by Caracci- 
olo et al. in a nice recent work [ p3| with our alternative 
theory to verify if they lead to better data collapse than 
when used to test the standard DDS equation (hopefully 
without having to introduce strong corrections to scaling 
and providing good order-parameter scaling) . 
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